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Abstract 

We discuss gauge freedom within the scope of the enlarged Helmholtz theorem 
and Neumann-Debye decomposition and then demonstrate its realization for 
the multipole expansion of a electromagnetic current with distinguished toroid 
moment family. The exact solution to the latter problem was obtained in 1974, 
but answers to some purely mathematical questions raised by LB. French and 
Y. Shimamoto Q about 40 years ago are given in this paper for the first time. 

I. INTRODUCTION 

In considering problems of mathematical physics with definite spatial and/or dynam- 
ical symmetries, one commonly uses various decompositions of vector fields over scalar po- 
tentials. These decompositions supplement the famous Helmholtz theorem and reduce its 
"gauge freedoms" . 

Let 3 V : r G R 3 — > V(r) e R 3 with all "good" properties. We may represent the given 
field V in different ways 

V -> {14, V y , V z } -> {<p, A} divA=Q -> {<p, if), x}- 

In any case, it is necessary to constrain superfluous components, if they take place, by 
introducing conditions similar to divA = 0. 

The last variant, diffeomorphic scalarization of a vector field, is the most economical and 
convenient approach to vector boundary-value problems of mathematical physics. But to 



use this approach, we have to be able to invert decomposition formulas, i.e. to deduce the 
integral representations of scalar potentials through the original vector field. 

Approach. To obtain a vector function from a (pseudo) scalar function set (p(r),ip(r) 
and x( r )) one must act on them by some vector operator F(r, V). Besides r and V them- 
selves, we may from them construct three simplest operators: L := — r x V, N : = 
V x L and M := — r x L. One can see their symplectic nature because of their correspon- 
dence to the frames of reference in a phase space: r, k, r x k, r x (r x k) and k x (k x r). 
The trio of vectors k, L, N and r, L, N, being immersed into the spaces i?| and Rp, 
respectively, form orthogonal bases in them, which are important for different applications. 

One may verify the following projection and commutation properties of F: 

[L,r 2 } = [L,p 2 ] = 

[L 2 , r] — i[r x L] - i[L x r] 

[L 2 ,k] = i[k x L] -i[L x k] 

rL = VL = rN = LN = LM = ML = 0, curlN = -LA; 

[L, A] = 0, [N, A] = 0, [M, A] = -6V, etc. 

[u, Vfc] = — 5 ik , [r-i, L k ] = —eikjfj, [Vj, = — ej fc jVj, 

[Vi, N k ] = VjVfc - A5 ik , [n, M k ] = nr k - r 2 5 ik 

[Li, M k ] = e ikj rj - r 2 e ikj Vj. 

It is taken into account that 

N = -rA + V(rV) + V, 
r x N = -(1 + rV)L = -L(l + rV), 
curlN = -LA, etc. 

Note that in the space Rj>, L and N form the following algebra of differential operators 
[Li, Lj] = 6ij k L k , [Li, Nj] = €ij k N k , [Ni, Nj] = —€ij k L k A. 
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After the rescaling of N — > N/\/S over a supporting function space, we acquire the 0(3, 1) 
Lie algebra representation in terms of the F(r, V) operators in R%.. They act in the 0(3) x 
0(3) function space suitable for arranging the multipole phenomenology in electromagnetic 
theory (see e.g. 0, |J, ||). The operators L and M obey the same algebra in the space of 
wave vectors that is Fourier-conjugate to -RjL In the preceding paper J7J the emphasis 
has been made on the inversion of different formulas for decomposition of vector fields in 
the mathematical aspect. 

In the physical aspect, two vector-potential formulations of the theory of continuous 
media with taking into account both magnetic and electric toroid polarizations |§, (see also 
||, part II) are published for the first time. Here we consider the mathematical underlying 
reason of uniqueness of division of the transverse electric distribution density Ei m (k 2 , t) into 
two independent multipole specimen 0: 

E lm (k 2 , t) = Q lm (0, t) + k 2 T, m (k 2 , t), (1) 

where Qi m (0, t) are the time-derivative of the Coulomb (charge) multipole /—moment and 
T; m (k 2 , t) are the toroid multipole form factor of /—order. Note that a secularized relation 
when neglecting toroid contributions and known as the Siegert theorem may be correct 
for low-symmetric electromagnetic systems only. The mathematical question of condition 
responsible for the possibility of identification of multipole moments (the leading ones for 
each given /) in the transverse and longitudinal parts of current raised in |J will be tried 
to clear in this paper. From our representation theorem (see also remark II) it follows that 
definitions of multipole moments are unique. 

Recall that the procedure of multipole decomposition of the current and field densities in 
electrodynamics in fact corresponds to the description of properties of a system by a set of 
numerical characteristics which are assigned to a point, "center" of the density distribution 
of the system considered. In this case, for the poloidal and potential parts of the current 
an additional connection arises between its longitudinal and transverse components due 
to degeneracy of boundary conditions of longitudinality and transversality in r-space at the 
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self-similar shrinkage of the definition domain of current to the chosen center. To prove this 
statement, we could probably use the transfer technique of boundary conditions ||], @. It 
might have lead to the separation of multipole moments called in H the toroid ones *. But 
ways of this kind is very difficult in the general framework of the distribution theory Here 
we use the simple concrete approach. 

The main feature of the multipole expansion procedure is a special choice of basis func- 
tions, which in actual practice ensues the rapid convergence of multipole series. This circum- 
stance forces us to weaken requirements of the usual Helmholtz theorem and, respectively, 
to take account of a gauge freedom extension. In sections 1 and 2, we discuss gauge free- 
doms in the Helmholtz and Neumann-Debye decomposition. In section 3, we turn to their 
realizations within the multipole expansion of the electromagnetic current. 

II. THE ENLARGED HELMHOLTZ THEOREM 

We begin our consideration with the Helmholtz Decomposition: 
< VV with properties of single-valuedness, continuity, boundedness or convergence 

. , k 

V < , e > at r -> oo 

1 ' /pZ-\-£ 



in the space may be represented in the form [see e.g. ||12 |] 



V = V(f + curlA with divA = 0. 

Indeed, the theorem requirements are proved to be sufficient in order to reexpress tautolog- 
ically the given vector field V through its divergence and vorticity. The explicit realization 
of the theorem could be attained due to the following operations 



divV 11 = A(p, <p = A _1 divV 
where A -1 := — 



d 3 r' 



'qc_r3/{o} 47r|r — r'| 



*In fact, their exact title should be the poloidal moments, see ||, 
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curlV" 1 = curl curlA = -AA, A = -A ^urlV. 

So we have 

V = V 11 + V 1 = VA -1 divV - curlA _1 curlV. > (2) 

Form ([!]) often produces misunderstanding (e.g. that the representation of divV and 
curlV is equivalent to the representation of V itself. We discuss here to what extent single- 
valued is this representation. In fact (1) contains the evident "gauge freedom": 

(p = (p + u, (Alu = 0), A = A + Vw. 

If we remove the demands of topological triviality Q or/and the boundedness of functions u> 
and w, then as an example important for physical application, we may represent the gauge 
freedoms in the form of special additional functions to V: 



0^V N = curlL( i/ r j +1 )F, m EE( J 



-(I + i)Vr*y /m N 
Vr-'-iy, m , 



which have a nonzero finite value all over the space except r — > oo and r — > 0, respectively, 
as far as 

divV N = curlV N ee 0, in R 3 /S*_> 00 and R 3 /{0}. 

Therefore these functions [] cannot be represented by the usual Helmholtz decomposition, 
and manifest the gauge freedom of its enlarged formulation. 

Remark I. Note that functions V^r are longitudinal and transverse simultaneously, since they 
represent the vector solutions of the Laplace equation (A = graddiv — curl curl). Moreover, in 
our context it is important to emphsize that they are topologically equivalent to the poloidal 



^Note that similar functions may be generated using the common relation 

curl(r x V)r K Y lm = -(« + l)Vr K Y lm + (« - + I + l)rr K - 2 Y lm , 
and taking into account that 

k = I and k = — I — 1 



(see [14 1, Appendix A). 



(meridional) harmonics on the toruslike surface covering the whole space of R 3 except one deleted 
axis. 

Thus, under this gauge freedom the Helmholtz decomposition takes the following alter- 
native forms 

V = V 11 + V ± = 

= V ( A-MivV - (/ + 1) YX°iry + CLr^Yin) + curlA^curlV = (3) 

V Im ) 

= VA-MivV + curl ( A-WV + L"£[C lm r l + C| m r- , - 1 ]^ m ) . 

V Im ) 
III. THE NEUMANN-DEBYE DECOMPOSITION 



The well-known mathematical physicist W.M.Elsasser in [15] has already observed 
that every vector field of the form VxVxr^ + Vx rip, where ip and x are anv scalars, is 
solenoidal. In paper |TH] it has been shown that if divV = in R 3 , then for every choice of 
the origin there exist unique scalars ip and x such that V = hip + N% while ip and x average 
to zero on every spherical surface concentric with the origin. The complete theorem of the 
possibility of decomposition of V(r) in terms of scalar functions reads as follows: 

Representation Theorem. Given a region £1 C i? 3 \{0}, with a regular boundary 
and R 3 — vector field, V : r G Q — > V(r) G R 3 . Then, there exist three scalar functions 
ip(r),ip(r) and x(r) on Q which define this V [|l7j]. ■ 

The most used decomposition, which we call the Neumann-Debye one, has the form 

< V(r) := Vip(-) + curlr'0(-) + curl curlrx(-) = Vy9 + Lip + N%. (4) 

Here ip and x are the so-called Debye potentials and tp is the usual (electric) scalar one. We 
found fundamental solutions of the inversion problem of (1) in the form 

divV = -> tp = A _1 divV, 

LV = L 2 ip -> ip = -L- 2 LV = L" 2 r curlV, (5) 
rV = (rV)cp + L 2 X -> X = L~ 2 (rV) A"MivV - L" 2 (rV), 
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where |T8 

dw' 



L 2 



f dw' , A 
/ — — mil — r ■ r ). > 



Remark II. Eigenfunctions of the square of the operator angular momentum iL : (iL) 2 = — L 2 
are usual spherical functions satisfying the equation 

L 2 Y lm (f) = -1(1 + l)Y lm (f) 

The corresponding Green function for this equation can be found with the help of the known 
Mercer theorem ( flSfl, v.l) which in the given case yields 



Remark III. It is well-known that the gauge freedom of (||) is the following: 
Requiring ip to vanish on the boundary, and ip, x n °t to contain spherically symmetric components 



dw = dw x = 0; 
s 2 Js 2 

we put these functions in one-to-one correspondence to V. ■ 

Uniqueness Theorem. If a vector field V (with the properties determined in the 
theorem (I) ) is defined on every S 2 in some range r < r < r\ and in that range V r = while 
Vg(r,9,ip) and V^(r, 0,(p) are bounded for each fixed r and are continuously differentiable 
except possibly at 9 = and 6 = 71 and further divV = and cur IV = 0, then V = 0. 



Our inversion formulas (3) demonstrate that immediately (cp. with [16|, p. 383, where 
the condition LV = has been used instead of our curlV = 0). ■ 

Now we have to reconstruct the representation of x such that it depends on curlV 
and divV only. Really, the latter quantities have the physical meaning but not the radial 
component of V. Moreover, we may expect \ not to depend on divV generally because this 
potential defines the transverse part of the vector field V. Nevertheless, because of gauge 
freedom, the situation is not so simple as it seems to be. 

Indeed, we may substitute the Helmholtz decomposition (the last expression in (^)) into 
the term with rV and see that x takes the form 



X = L^LA^cmlVj^L-^C^Vr 1 + C^Vr"'- 1 ]!^. (6) 



Irn 



Further our vector field corresponding to the gauge freedom may be transformed as 

curlLL- 2 (rV)r^ m = (I + l)curlr'LL- 2 ^ m = 
= (l + l)cur\r l LL~ 2 j^-^Y lm = jcurlLr^ m = --j^Vr l Y lm . (7) 

So, we found fundamental solutions of the inversion problem of @ in the form 

divV = A(p -> (p = A -1 divV, 

LV = L 2 i) -> V = -^ _2 LV = L~ 2 r curlV, (8) 
rV=(rVV + L 2 x -> X = r 2 (rV) A _1 divV - L~ 2 (rV), 

= L^LA^curlV 

Irn 

with taking into account the gauge freedom in (§). ■ 

IV. THE MULTIPOLE EXPANSION OF THE 
ELECTROMAGNETIC CURRENT 

Now we compare our abstract exercises and the procedure of multipole expansion of 
electromagnetic current on the basis of the vector Helmholtz equation solutions constructed 
through the Neumann-Debye decomposition (|8p. The multipole representation of current 
J(r) by (H) may be obtained by the standard expansion of three scalar densities into series: 



~ ^2ji{kr)Y lm Qi m (k 2 ,t); 

Irak 

Tp~J2ii(kr)Y lm M lm (k 2 ,t); (9) 

Irak 

X~J2ji(kr)Y lm E lm (k 2 ,t). 

Irak 

Therefore the multipole representations of the transverse part of current (divJ = 0) are 
determined by the magnetic form factors Mi m (k 2 , t), and the transverse electric contributions 
Ei m (k 2 , t), and the scalar part of current (curlJ = 0) are expressed due to the 4— current 
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conservation law (div J = —p) through the Coulomb (charge) multipole moments Qi m (0t) 
and their mean 2n— power radii 

oo yjln 



Qim(k 2 , t) = Q lm (0, *) + £ TT^-^ ' *)■ 
As a result, J may be represented as || 



n=l n[ 



J(r, t)= (2tt)- 3 E / dk(-ik) 1 V ; / {-»fc Lf^Y^M^k 2 , t) 

-ikcniYLf^YU^EUk 2 , t) + I VMkr)Y lm (r)Q lm (k 2 , t)}. (10) 

We repeat here the procedure used firstly in || for the ascertainment of the exact struc- 
ture of transverse electric contributions Ei m (k 2 , t). To this end, it one should rewrite ex- 
plicitly our basis functions in terms of the vector harmonics 

curlL/,(fcr)Y, m (r) = (21 + l)^ 2 {f^h^y/T+lY U -i m {r) + 

+fi + i(kr)VlY ll+lm (i)Y, 

V/,(Ar)Y, ro (r) = (2/ + ^-^{^(A^Y^r) - 



-fi +1 ( kr ) VTTlY a+lm (f) } , 

where Y^_i m (r) is a harmonic polynomial function defined as 

r'- 1 Y /z _ lm (f) = 1 V r'F im . 
^(2/ + 1) 

As is obvious, in the wave-length approximation, the leading contributions in the latter 
expressions are delivered by the vector harmonic functions fi-\Yu-i m from which it follows 
that 



curlL/,(£r)}Ur) « fc ^, <J(l + l)/JV/,(Ax)Y, m (r) 



j-Vr^ m (r) + 0[(k) m ] (11) 



(2/+ 1)!! ' 

It is the relation that permits us to identify the leading term in Ei m (k 2 , t) with time- 
derivatives of Qim(0, t). However all functually independent contributions in Ei m {k 2 , t) give 
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the so-called toroid moments and their 2n— power radii 

1(1 + 1) 1 

J(r,t)d 3 r. 



Tt\0, t) = -^L Jr^ 



V*/(* + 1) 

21 + 3 



Y^_i m (r) H — Y n+lm ( 



As the toroid moments have a distinct geometrical meaning (diverse details and represen- 
tations are given in ||, [f|, |§, [K| and see also ]ll[), the rejection of TJ m (i) is generally 



invalid like it was done in the Siegert theorem Ei m (k 2 , t) — >fc_,.o Qim(®, t). Neglect of TJ m (t) 
in comparison with Qi m {t) is analogous to the neglect of a higher multipole moment (con- 
tributions of highest symmetries of a given system) in comparison with the lower ones, 
which is of course permissible only when the lower moments (low symmetries) of this sys- 
tem do exist. So, the strict theorem determining the electric part structure has the form 
(|l|) Ei m (k, t) = Qim(0, t) + k 2 T/ m (k, t), and its validity and uniqueness rely on the gauge 
freedom which has been obtained for the enlarged Helmholtz theorem and transferred to 
Debye potentials (compare (|8|) and (|TTD). 

Moreover, by using the exact relation (|T[), it is not so hard to find an expression for the 
complete parametrization of the current in terms of generalized functions |§ 



oo m=l oo fcyj _i_ 1 \ 1 1 / /| 

^ m tl^ 2"n!/(2Z + 2n + l)!!V2Z + l 1 lm wv ' KJ 



+ [Q lm (t)S nfi A- 1 - T£ n) (t)]curl(r x V)A n 6 lm (r) - I Q^\t) VA^Mr)}. 



(12) 



Remark IV. However, the expansion via of spherical harmonics emerges rapidly convergent. 
We remind (see e.g. (2j, p. 806) that if the field decreases for large distances very slowly, slower 
than r -2 , the divergence and curl of the vector field considered are assigned arbitrary independent 
values. Conversely, if we know that J vanishes identically outside some source radius R, VJ and 
V x J are no longer independent of each other. As far as the expansion via spherical harmonics is 
rapidly convergent, it is realized in the latter representation (|l^) immediately. 



^The contribution of I = 0, n = to the last term is forbidden by total charge conservation, 
whereas other terms contain no contribution of I = 0, formally, owing to (r x V)<5(r) = 0. 
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Thus, we strictly demonstrated that the gauge freedom in division of the electromag- 
netic current into the transverse and longitudinal parts leads to the fact that the multipole 
contributions to the transverse part of current E[ m (k 2 , t) are represented in the form (|IJ) and 
its leading terms may be identified with Qi m (0t) from the longitudinal part of current for 
all I. Note that, since the coefficient C\ does not depend on the wave number k, we can use 
hereafter the Lorentz gauge condition in the calculation of the vector potential. 

V. CONCLUSION 

The representation of \ in the Neumann-Debye scalarization already assumes that the 
prohibition of the electric type of radiation imposes some conditions both on curlJ and divJ. 
But we could not reveal the ones due to their non-division in the scalarization mentioned. 
Exploitation of the enlarged Helmholtz theorem for this operation has inserted the extended 
gauge freedom (let us recall that the Neumann-Debye representation the gauge freedom 
reduce to functions of the scalar argument |r| only). It is just this freedom, consideration 
of which made it possible to identify the coefficients of leading order of the expansion of 
transverse and longitudinal electric parts of the current! 

The form of the expression ([!]) shows the possibility of compensation of the electric type 
radiation if the toroid and charge moments are switched on as "anti-phase" ones || (see 
also H). 
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